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On the Expectations of Maxima of Sets of 
Independent Random Variables 
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Abstract 

, ^ ■ Let X 1 , . . . , X k and Y 1 , . . . , Y m be jointly independent copies of random 

variables X and Y, respectively. For a fixed total number n of random vari- 
ables, we aim at maximising M(k, m) := EmaxjX 1 , . . . , X k , Y l , . . . , Y m } in 
k = n — m > 0, which corresponds to maximising the expected lifetime of 
an n-component parallel system whose components can be chosen from two 
p> | different types. We show that the lattice {M(k,m): k,m > 0} is concave, 

give sufficient conditions on X and Y for M(n, 0) to be always or ultimately 
maximal and derive a bound on the number of sign changes in the sequence 
M(n, 0) — M(0,n), n > 1. The results are applied to a mixed population of 
Bienayme-Galton- Watson processes, with the objective to derive the optimal 
initial composition to maximise the expected time to extinction. 

Key words and phrases. Parallel systems, expected lifetime, branching processes. 
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\6 
o 

o\ '. 1 Introduction 

O 

Consider a simple ra-component parallel system (called an n-assembly in what fol- 
^ ■ lows), whose component lifetimes are independent non-negative random variables 

Xi, . . . , X n with finite means. The n-assembly fails when all the n components fail, 
and so its expected lifetime is given by Emax{X 1; X 2 , . . . , X n }, which is a standard 
measure of the system's reliability. 

Such quantities and related characteristics of complex systems often appear in 
literature on reliability. For instance, comparisons and bounds for expected life- 
times of series and parallel systems are considered in [S], while assemblies whose 
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components' lifetime distributions can depend on the amount of "investment" are 
discussed in [7]. Expected lifetimes conditional on survival up until time t are in- 
vestigated in pQ, while assemblies of dependent components with interchangeable 
lifetime distributions are analysed in [9] . More recently, upper and lower bounds for 
the expected lifetimes of n-assemblies of components with (possibly) non-identically 
distributed lifetimes in terms of the expected lifetimes of assemblies with identically 
distributed components' lifetimes were given in [S]. 

In the present study, we assume that, when building an n-assembly, one may 
choose components from d different types that have different lifetime distributions. 
Should one choose all n components of the same type (unmixed assembly) or should 
one mix? In general, the answer will depend not only on the distributions of the 
random lifetimes of the components, but also on n. In this paper we shall examine 
these relationships. 

Questions of such kind arise naturally in reliability problems such as server fail- 
ures, where information is duplicated on multiple hard disks. In conservation biology, 
the "Single Large or Several Small" problem provides another application. Here one 
should decide between the creation of a single large protected habitat or several small 
ones. This may be modelled by n-assemblies since the lifetimes of individuals in the 
species in different small habitats may be modelled by random variables Xi with 
distinct (due to differing environmental conditions) distributions. In contrast, a sin- 
gle large habitat with the same lifetime distribution for all individuals in the species 
corresponds to an n-assembly composed of identically distributed random variables. 
Furthermore, the question of whether mutation is inherently advantageous because 
it leads to diversity may also be analysed in a similar way. 

In this paper, we first show that the expected lifetime M(k\, . . . ,kd) of an as- 
sembly consisting of ki components of type i = 1, . . . , d, is a concave function of 
(ki, . . . , kd)- This implies that one can easily find the optimal n-assembly in the 
general case using the discrete steepest ascent method in linear time in dn. After 
that we focus on the case of two component types, with a view to maximising the 
life expectancy of the system, noting that generalisations to more than two types are 
straightforward and are omitted for the sake of simplicity. We denote the two ran- 
dom lifetimes by X and Y and their respective distribution functions by F and G, 
always assuming that ~EX + KY < oo, and assume that X\, X 2 , . . . are i.i.d. like X, 
and Y"i, Y2, . . . are i.i.d. like Y, with all the X^s and Y/s being jointly independent. 
Our main object of study is 



We classify the possible behaviour of these quantities, e.g. when it is advantageous 
to always choose one type and when it is the right choice eventually, giving sufficient 
conditions for these. We also give a simple bound for the number of sign changes of 
M(n, 0) — M(0, n) and apply the aforementioned results to the case of Bienayme- 



M(k, m) := Emax{Xi, . . . , X k , Y u . . . , Y m } 
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Galton- Watson processes, with a view to selecting an optimal initial composition in 
order to maximise the expected time to extinction of the entire population. But first 
we present a couple of simple examples that show some possible types of behaviour 
that can be displayed by M(k,m). 

Example 1. Suppose we can choose from two types of components: one with life- 
times X ~ C/[0, 1] (i.e. distributed uniformly on [0,1]) and the other with deter- 
ministic lifetimes Y = e = const. Clearly, if e > 1 then one should always choose 
components of type Y . In general, we will refer to a situation where the optimal 
(with the longest expected lifetime) n-assembly is always composed of one partic- 
ular type of component as the dominant case, and to that component type as the 
dominant type. 

Now assume that < e < 1. The choice is between selecting all n components 
of the first type or n — 1 components of the first type and one component of the 
second type (there is no point in taking more than one type two component to 
extend the expected lifetime as all such components will fail at precisely the same 
time). Integrating to obtain the corresponding expectations gives 

r 1 1 

M(n,0) = y (l-x n )dx = 1-— — p 

f 1 1 £ n 

M(n - 1, 1) = e + / (l - a;" -1 ) dx = l-- + —. 

J £ ' n n 

We conclude that the assembly corresponding to (n, 0) yields a higher expected 
lifetime than that for (n — 1,1) for n such that (n + l) l l n < e~ l . If e < 1/2, 
this will always be true, and then we have a dominant case. Since e < 1 and 
(n + l) 1//n — > 1 as n — > oo, we see that, regardless of e, M(n, 0) > M{n — 1,1) 
holds for all large enough n. In general, if unmixed assemblies {X\, . . . ,X n } are 
advantageous for all large enough n, we say that X is ultimately dominant. We 
will see (in Theorem E]) that this behaviour takes place whenever one distribution 
function ultimately dominates the other. 

Since (n + l) 1 /™ is a decreasing function of n > 1, we see that the values of n for 
which [n + l) 1 /" > e _1 , i.e. for which {n — 1, 1) is preferable to (n, 0), form either an 
empty set or a finite sequence of successive integers. In the latter situation, since 
the optimal n- assemblies are composed of more than one component type, we say 
that we have a non-dominant case with the number of dominance changes (i.e. sign 
changes of M(n, 0) — M(0, n)) equal to 1. 

We will see that this number admits a simple upper bound in terms of the 
distribution functions of X and Y (Theorem [7]), but can actually be infinite, as the 
following example suggested by A. Sudbury shows. 
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Example 2. Let X and Y have distributions F(X = 2i) = -pi and F(Y = 2i+l) = q t , 

1 = 0, 1, 2, . . . , such that 

Pj ■= Pj + Pj+i + ■■■ = Qj+i + Qj+2 + ■■■ = 2 _2J , j > 2. 

Set A(k, m) := max{Xi, . . . , X k , Y\, . . . , F m }, so that M(fc, m) = EA(k, m). 

We will show that there exist two increasing infinite sequences {hi} and {/«} such 
that M(ki,0) > M(0,ki) and M(h,0) < M(0,k) for all sufficiently large i. 

Indeed, taking ki := \2 2 'i~ 2 ], where \x] = min{A; G Z : k > x}, we see that 

F(A(k h 0) < 2i) = (I — i*)* 4 = (I - 2- 21 ) 221 ^ 2 = exp{-(l + o(I))r 2 } - 1, 
P^, 0) < 2{i - I)) = (I — 2 - 21 ") 221 ^ 2 = exp{-(l + o(I))2 2l ~V 2 } -> 

as i — > oo, so that M(/cj,0) ~ 2(i — f) (here ~ means asymptotic equivalence: 
u n ~ as n — > oo iff u n /v n — » I), whereas 

P(^(0, fci) > 2i + 1) = 1 - (1 - Pi-i)* 4 

= I - (I - 2 - 21 ") 221 ^ 2 = I - exp{-(l + o«)2 2l ~V 2 } - I. 

Thus M(jfei, 0) < M(0, fcj) for all sufficiently large %. 

Now let U := 2 2 (i + iy 2 . Then it can be shown in a similar fashion that 
M(/j, 0) > M(0, /j) for all sufficiently large i, so that we have infinitely many domi- 
nance changes in this example. 

2 Concavity of the expected lifetimes and 
dominance 

When maximizing the function M(k\, . . . , kd), it is helpful to know that it is concave. 
In particular, this guarantees that, moving in the direction of the steepest ascent 
one can find the optimal composition in linear time in dn. 

Theorem 3. Let Xi,...,X d be random lifetimes with corresponding distribution 
functions F±, . . . , Fa. Then the function 

POO 

M(x) := / (I - F^ ■ ■ ■ F d (s) x «) ds, x = (x u ...,x d )e R d + , 
Jo 

is concave. 

Clearly, the function M{x) coincides with the corresponding systems' expected 
lifetimes at the integer lattice points. 
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Proof. We will begin with the obvious observation that, for any x G Mi, 

d 

M(x) < M(\ Xi ],. . . , \x d ]) < Y^\xjWCj < oo. 

Fix arbitrary x, y E M.+ and let F(s) := F^(s) ■ ■ ■ F* d (s), G(s) := Ff (a) ■ ■ ■ F y d d (s). 
It suffices to show that the univariate function 

poo 

M 1 (z) := M(zx + (1 -z)y) = / (1 - F(s) z G(sy- z ) ds, ze [0,1], 

Jo 

is concave. To justify differentiation inside the integral, we first fix an arbitrarily 
small 5 > and show that the integral of the integrand's derivative converges 
absolutely and uniformly in z G (5,1 — 5). We have 

poo a poo 

/ — (1 - FisYGis) 1 -') ds= / \\nF(s) - \nG(s)F(s) z G(s) 1 - z \ds, 
Jo dz y Jo 

where the integrand is dominated by the function 

\ In F {s)\F {s) s + \ In G(s)\G(s) 5 

(using the natural convention • InO = 0), which is clearly bounder on (0, oo) and 
so is integrable there since | \nF(s)\F(s) s ~ 1 — F(s) as s — ► oo. 

In a similar way one can verify that Mi is twice differentiable and 

roo o2 

M»(z) = J o —(l-F( S yG( S y- z (s))ds 

POO 

= - (lnF(s) -\nG(s)) 2 F(s) z G( S y- z ds < 
Jo 

unless F(s) = This shows that Mi and hence M are concave. □ 

Now we turn our attention to the case when one type is always dominant, i.e. 
one has M{n, 0) > M(k, n — k) for any n G N, < k < n. 
It is obvious that if 



F(s) < G(s) for all s G 
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then we have a dominant case. One can weaken (pTJ) by requiring that one distribution 
always has a heavier integrated tail than the other. 

Theorem 4. (Theorem [36 of [3J) Let X and Y be random lifetimes with df's F 
and G respectively. If 

/oo 
(F(t) - G(t)) dt<0 for all s G R+, (2) 

then M(n, 0) > M(k, n - k) for all n G N and k G [0, n) . 
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Recall that ([2D defines the so-called increasing convex ordering on the set of 
distributions on K + (see e.g. p. 11 in [1 lj ) , whereas ([I]) is the standard stochastic 
ordering. 

The next example shows that the above condition can be necessary in some 
special cases. 

Example 5. Let X be uniformly distributed on [0, 1] and Y uniformly distributed 
on [a, a + e] with 0<a<a + e<l. Then clearly 



«i 

M(n, 0) - M(0, n) = / (l-x n )dx- 



o 



(/ / ( 1 (- ) ) (If 



a 



= ( 1 - £ »( 1 -^tt)- 

Hence M(n, 0) > M(0,n) iff a/(l - e) < 1 - (n+ l) -1 . This is the case for all n 
(and then the first type is dominant) iff a/(l — s) < 1/2, and so this condition is 
necessary for the first type to dominate. On the other hand, by Theorem H] for the X 
to be dominant it suffices that (j2J) holds, which is equivalent to the condition that the 
area of the upper of the two triangles formed by the graphs of the functions F and 
G is no less than the area of the lower one. Since the two triangles are similar, that 
means that the abscissa of the point of the intersection of the two graphs does not 
exceed 1/2, i.e. a/(l — e) < 1/2. Thus we see that condition (j2J) which is sufficient 
for dominance of X is equivalent to the condition that was shown to be necessary 
for the dominance. 

There is also a simple sufficient condition which ensures that one type will even- 
tually be dominant. 

Theorem 6. Let X and Y be random lifetimes with the corresponding distribution 
functions F and G, and let there exist an s £ ^+ such that F(s) > G(s) for all 
s > s . Then there exists an n e N such that, for all n > n and 1 < k < n, 
M(k,n — k) < M(0,n), i.e. Y is ultimately dominant. 

Proof. We will prove the result in the continuous case, an extension to the general 
case being obvious. 

Let s = inf{s : F(t) > G(t), t > s} and F(s ) = G(s ) =: e~ a . Then 

POD 

M(k,n-k)-M(n,0)= / (F(s) k G(s) n - k - G(s) n ) ds 



o 



F(s) x 



o 



o Js 



6 



Clearly, |Ji| < s e an . Now set H(s) := lnF(s) — lnG(s) > for s > s and choose 
si > so such that G(si) =: e~ b > e~ a . Then 

poo poo poo 

I 2 = G(s) n (e kH{s) - 1) ds > / G(s) n (e kH{s) - l) > fee" 6 ™ / ds. 

./so J Si J Si 

/oo 
H(s) ds > Soe~( a ~ b ^ n for all large enough n, which 

clearly implies that, for all k > 1, n> n , 

POO 

I 2 > ke~ bn / H(s) ds > s e- an > \h\, 

J Sl 

as required. □ 

Using Theorem [6], one can easily find examples of non-dominant cases. It suffices 
to have EX > EY and F(s) > G(s) for s > s E R+, so that M(l, 0) - M(0, 1) > 0, 
but M(n*,0) — M(0,n*) < for some n* G N (in other words, at least one sign 
change of M(n, 0) — M(0, n) occur between n = 1 and n = n*). In such situations, 
there typically exist n and k* e (0, n) such that M(k*,n — k*) > M(k, n — k) for all 
< k < n, i.e. the mixed assembly (k*,n — k*) is optimal. 



3 The number of sign changes of M(n, 0) — M(0, n) 

One can obtain a useful upper bound for the number of sign changes in the sequence 
M(n, 0) — M(0, n), n > 1, by generalizing the Descartes rule of signs. Clearly, this 
number does not exceed the number of zeros of the function 

poo poo 

ip{x):= (1-F(s) x )ds- (l — G(s) x ) ds, x > 0. 
Jo Jo 

In order to find a bound for this number, we will employ a generalisation of the 
Descartes rule of signs [5] . 

One says that an ordered system {ui,u 2 , ■ ■ ■} of functions defined on / C R 
satisfies the Descartes rule of signs on / if the number of zeros (with multiplicities) 
of a linear combination Y2 a j u j( x ) °f t nese functions is less than or equal to the 
number of variations of strict sign in the sequence {a\, a 2 , ■ ■ ■} of the coefficients in 
the combination. 

A system of functions {u\, . . . ,u n } defined on I C K. is said to be strictly sign 
regular if, for any fixed k < n and any X\ < x 2 < ■ ■ ■ < x m , m > k, with all 
Xi El, ? < m, all minors of order k of the matrix {uAxi)).^ have the same 

strict sign. For example, systems of the form {x^ 1 , . . . , x 13 "} and {e" 1311 , . . . , e _/3n:r } 
with < Pi < p 2 < ■ ■ ■ < p n are strictly sign regular on M (see e.g. p. 34 of [5]). 
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Proposition 2.6 of [5] asserts that a system of functions {ui, . . . ,u n } satisfies the 
Descartes rule of signs iff it is strictly sign regular. 

Given two distribution functions F and G on M + , let 

/2(A) := / [F- 1 ^-*) - G-\e- t )] dt 

J A 

be a signed measure on R + . Let S 1 " 1 " be a positive set and S - = M + \ 5+ a neg- 
ative set from the Hahn- Jordan decomposition \x = li + — lT . We assume (for the 
moment) that S + has only finitely many connected components , . . . , and 
that, without loss of generality, G S^. By S^, . . . , 5*~_ we similarly denote the 
connected components of S~, so that x± < x\ < x 2 < x 2 < ■ ■ ■ for any xf G Sf 
(note that < v + — i/_ < 1). Let v = v~ + v + and introduce sets Qi, . . . , Q v by 
putting Qi := S*, Q 2 '■= S^, Q3 := S 2 , Q4 := S 2 etc. The quantity v may be 
called the number of intersections of F and G since v + and v~ are none other than 
the numbers of connected components of {s : F(s) < G(s)} and {s : F(s) > G(s)}, 
respectively. When S + has infinite number of connected components, we set v = 00. 
This was precisely the case in Example [5] from the Introduction. 

Theorem 7. Let X and Y be random lifetimes, F and G being their corresponding 
distribution functions. Then the number of sign changes of M(n, 0) — M(0,n), 
n > 1, does not exceed v — 1. 

Proof Noting that X = F~ l (U) andF = G'^U), where U ~ U[0, 1] and H-\u) : = 
inf{t : H(t) > u} denotes the generalized inverse of the distribution function H, we 
can write ^ 

V»(x)= / F-\s 1/x )ds- I G-\s 1/x )ds. 
Jo Jo 

Letting s 1 ^ = e"*, we see that ip(s) = xipi(x), where 

^i(rr) = / e- xt fi(dt) = V / e~ xt fi(dt) 
Jo j=1 Jqj 

= ^(-i)^ 1 / e ~ xi IM^)I = $^ a i Ee_a:BJ 

i=i ^ i=i 

with ctj = fi(Qj) (note that sign (a j) = (— l)- 7 ^ 1 ) and Bj = H~ (Uj), where 

Hj(t) :=- f\ Qj (s)^ds) 
a j Jo 

are distribution functions supported by the corresponding sets Qj, j = l,...,u, and 
Ia denotes the indicator of A, while Uj ~ U[0, 1] are independent random variables. 
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To show that {uj(x)} := {Ee _irSj , j = 1, . . . , i>} is a strictly sign regular system 
of functions on R + , first note that {vj(x)} := {e~^ jX } is such a system provided that 
< Pi < fa < ■ ■ ■ < P u (see p. 34 and Corollary 3.9 of [5]). From this it follows 
that, for a k < u, the sign of det (e j^-^ is one an d the same for any choice of 
< Pi < ■ • • < p k and < xi < ■ ■ ■ < x k - 

Now observe that a kth order minor of i<m admits an integral repre- 

sentation in terms of the corresponding minors of (uj(xj)) . <w . <m . Indeed, consider 
without loss of generality the determinant of (uj(xi)) t . <k for < %i < ■ ■ ■ < Xk- It 
has the form 



k 



sgn(a) J] Uj{x a{j) ) = II ^ X ^ B > = E ^ sgn(a) J] e~ x 

cr j = l a j = l cr j = l 



J ■ ■ ■ J ^sgn(a)(n^ CTW )^i(A)---^(A) 



where the summation is over all permutations cr of {1, ... , k}, sgn(a) is the signature 
of cr and we used the independence of Bj. To complete our argument, it remains to 
make use of the above remark concerning the sign invariance of the determinants of 
the matrix (e~P jXi ) and recall that, by construction, we have Pi < P2 < ■ ■ ■ < Pk in 
the integrand in the last displayed formula, due to Pj G Qj, j < k. 

Thus we have shown that our {uj(x)} is a strictly sign regular system and so, 
by the Descartes rule, ipi (and therefore ip) can have at most v — 1 zeros in (0, 00). 
Thus M(n, 0) — M(0, n) can have at most v — 1 sign changes, as required. □ 

Note that it follows from Theorem [6] that if the (possibly infinite) sequence of 
points of intersection of the graphs of the distribution functions F and G is bounded, 
the number of dominance changes will necessarily be finite. However, if this sequence 
is unbounded (so that v = 00), the number of dominance changes can be infinite as 
Example [2] in fact shows. 



4 Mean time to extinction of mixed branching 
populations 

In this section we consider an application of the above results to the optimal selection 
for populations of simple branching processes, where n-assemblies are composed of 
the times to extinction of the progenies of the individuals comprising the initial 
population. 

Suppose we wish to establish a colony populated by two species, Si and S2, with 
the aim of maximising the colony's expected lifetime. Given a fixed initial total 
number n of individuals in the colony, we can choose k individuals of species Si 
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and n — k individuals of species S^. Next we assume that time is discrete and that 
the reproduction of individuals in the colony is governed by independent subcritical 
Bienayme-Galton- Watson (BGW) processes specific to the species. 

Let £ > be an offspring random variable with a distribution pi = P(£ = i), 
i > 0, and £ nj j with n, i — 1, 2, . . . be independent copies of £. Further, for r G N, 
let Z< r )(0) := r and 

ZM(n-l) 

ZW(n):= ^ n>l. 

Then {Z"(n), n > 0} is a BGW process starting with r individuals (for r = 1, we 
will simply write {Z(n)}). 

Denote by := inf{n > 1 : Z^ r '(n) = 0} the time to extinction of the process 
{Z( r \n)} (we write T for T^). Observe that = maxjT 1 , . . . , T r }, where the 
T J 's are independent copies of T. Let f{9) := YHiLo V$ l De t ne probability generating 
function (pgf) of £ and f n (6) the nth functional iterate of /. Then 

P(T < n) = P(Z(n) = 0) = / n (0), n > 0, (3) 

and hence P(T( r ) < n) = f n (0) r and ET (r ^ = ^ =0 (l - / n (0) r ), where we put 
./n(n) 0. 

Now denote by {Zi(n)} a BGW process describing the growth of a population of 
type Si individuals (starting with a single individual at time 0), i = 1, 2, and by Ti 
the respective times to extinction. We will keep the notation /, f n for the offspring 
pgf and its iterates for species Si and denote the respective pgf 's for S 2 by g and 
g n . Then the expected time to the extinction of our colony initially consisting of k 
individuals of species Si and n — k individuals of species S2 is given by 

M(k, n — k)—E max{T?\ T^ k) } 

00 

= Emax{Tl, . . . ,T*,T 2 \ . . . ,7^} = ^{l- /,(0)V t (0) n - fc ), 

i=0 

where the T- are independent copies of Tj, i = 1, 2. The problem we want to consider 
is how to choose that maximises M(k, n — k) for a fixed n. First we obtain some 
preliminary facts. 

Theorem 8. For any two subcritical BGW processes {Zi(n)} and {^(n)} with 
offspring means /j,2 < fii < 1, t/iere exists an no such that /„(0) < (?„(0) /or a// 
n > no and therefore T\ is ultimately dominant, i.e. there exists an n\ such that 
M(n, 0) > M(k, n — k) for all n > n\ and 1 < k < n. 

Proof. From [TU] we know that, as n — > 00, 

l-/„(0)~/i?L!K) and 1-^(0)~^L 2 (^), 
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where Lj(s) are functions slowly varying at 0, i.e., for any c > 0, Lj(cs) ~ Lj(s) as 
s \ 0, i = 1,2. As is well-known (see e.g. Proposition 1.3.6(v) in [4J), if L(s) is 
slowly varying at then, for any e > 0, one has s £ < L(s) < s~ £ for all small enough 
s > 0. Therefore, for an arbitrarily small e, 



4 > (s lS2 y 



> I — e I —* °o as n — > oo. 



once si and s 2 are small enough. So, for < e < (ln/i 2 — ln/!i)/(ln/i 2 + ln/ii) we 
obtain 

1 - /n(0) »l M/fl) ^ ^l +£ 
l-<7„(0) ~ ^L 2 (/i£) V/i2" 

Hence / n (0) < g„(0) for all large enough n. Together with that means that 
-F(s) < (j(s) for all large enough s. Therefore 7\ will be ultimately dominant by 
virtue of Theorem El □ 

Note that, in the case of critical BGW processes with offspring pgf's / and g, 
a result analogous to the first part of the above holds true: if the variance of / is 
lower than that of g, then / n (0) > g n {0) for all large enough n (this follows, say, 
from relation (1) on p. 19 of |2J). 

The next results provides an upper bound for the number of dominance changes. 

Theorem 9. Let {Zi(n)} and {Z 2 (n)} be two subcritical BGW processes with re- 
spective offspring pgf's f and g and times to extinction T\ and T 2 . Then the number 
of times the difference f n (0) — g„(0), n = 1,2, ... , changes sign and the number of 
dominance changes for unmixed n-assemblies composed of copies of T\ and T 2 are 
no greater than the number of zeroes of f(6) — g(9), 9 e [0, 1). 

Proof. Suppose that f{9) > g{9) for 9 G [a, b) C [0, 1] and / m (0) > 5^(0) for some 
m such that / m (0) G [a, b). Then clearly / m+i+1 (0) > g m+ i+i{f)) for alH > such 
that / m+i (0) G [a, b). Now supposing that, for some c G (6,1], f{9) < g{9) for 
9 G [b, c) and / n (0) < g n (0) for some n such that f n (0) G [b, c), we likewise deduce 
that f n+ j + i(0) < g n+ j + i(0) for all j such that / n+J (0) G [b, c). It is easily seen 
that the above implies that, for any interval [a, c) C [0, 1] with at most one zero 
of f(9) — g{9), there will be at most one change of sign in the finite subsequence 
{/n(0)-</„(0):/ n (0)e[o,c), n G N}. 

Suppose that f{9) — g{9) has v* < oo zeroes in [0, 1). Then, subdividing [0, 1) 
into subintervals with one zero in each, it follows that the total number of changes 
of sign in the sequence {/ n (0) — fi'n(O) : n G N} will not exceed u*. 

Now let F and G be the respective distribution functions of T\ and T 2 , i.e. 
F(s) = /[ s j(0) and G(s) = <7| s j(0). Then the number of connected components of 
S + := {s : F(s) > G(s)} = {s : / w (0) > <? W (0),s G R + } and S~ := {s : F(s) < 
G(s)} = {s : /[ s j(0) < g[ s ](0),s G M + } will not exceed u*. Applying Theorem [7] 
completes the proof. □ 
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Suppose now that the dynamics of the populations of the two species are modeled 
by subcritical branching processes {Zi(n)} and {Z 2 (n)} with offspring pgf's / and 
g, respectively, such that 

f(9) — g(9) has at most one zero in [0, 1). (4) 

This is the case, for instance, when both pgf's are quadratic: f{9) = p + p\9 +P2O 2 
and g{9) = q + q x 9 + q 2 9 2 with p + p 1 + p 2 = q + q 1 + q 2 = 1. 

The following result describes the behaviour of M(k, n—k), < k < n, depending 
on agreement of the signs of the differences in respective offspring means and mean 
times to extinction. 

Corollary 10. Assume that subcritical pgf's f and g satisfy condition (jlj). 

(i) If {iii — /i 2 )(ETi — ET 2 ) > then we have a dominant case, with M{n, 0) largest 
when fii > fi 2 and M(0, n) largest otherwise. 

(ii) // (/ii — fi 2 )(ETi — ET 2 ) < then we have a non-dominant but ultimately 
dominant case with exactly one sign change in the sequence M{n, 0) — M(0, n), 
n = l,2,... 

Proof, (i) Without loss of generality assume that \i\ > \i 2 . Then, on the one hand, 
T\ is ultimately dominant by Theorem [HI so that 

/n(0) — <7n(0) < for all large enough n. 

On the other hand, due to condition (jlj), by Theorem [9] the sequence f n (0) — g n (0), 
n = 1,2,..., can change sign at most once. Therefore there exists an n Q > such 
that 

/n(0) - <? n (0) > 0, n < n , 
/n(0) - g n (0) < 0, n > n , 

which clearly implies that 

X}t/n(°)-<7«(°)) <0, 3>n , 

whereas for j < no one has 

W0) - 9n(0)) < J2 - 9n(0)) = ET 2 — ETi < 

n>j n>0 

by our assumption. Hence Ti is dominant by virtue of Theorem HI 

(ii) We again assume that \i\ > fi 2 , which implies that 7\ is ultimately dominant 
owing to Theorem [HJ However, 

ETi = M(l, 0) < M(0, 1) = ET 2 

by assumption, so there is at least one change of sign in M(n, 0) — M(0, n). Finally, 
as f(9) — g(9), 9 G [0, 1), has at most one zero, Theorem [9] asserts that there cannot 
be more than one change of sign in the sequence. □ 
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